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ABSTRACT 

The  concept  of  2D  elasticity  of  a graphene  sheet  together  with  the  idea  of  stiffness  of  a single 
sp3  bond  have  been  applied  to  theoretical  evaluating  elastic  properties  of  diverse  carbon  states. 
2D  elastic  moduli  have  been  extracted  from  data  on  elastic  moduli  of  crystalline  graphite. 
Stiffness  of  the  sp3  bond  has  been  estimated  from  data  on  the  elastic  modulus  of  diamond. 
Efficiency  of  Van-der-Waals  interaction  has  been  taken  from  the  elastic  modulus  C33  of 
crystalline  graphite.  Characteristics  of  single  fullerene  deformability  have  been  computed  by  the 
molecular  dynamics  method.  Theoretical  estimations  have  been  performed  for  single  molecular 
clusters,  pristine  fullerite,  HPHT  phases  of  polymerized  C60,  etc.  The  estimations  are  in  good 
agreement  with  experimental  data  on  elastic  properties  and  nanoscale  structure  of  carbon  states. 
The  approach  is  effective  for  establishing  interrelation  between  nanostructure  and  elastic 
properties,  for  prediction  and  classification  of  nanostructure  in  novel  carbon  materials. 


Three  types  of  hybridized  covalent  bonds  between  of  carbon  atoms  (sp3,  sp2  and  sp)  give  rise  a 
variety  of  macromolecular  clusters  and  carbon  nanostructural  solids,  which  are  essentially 
distinguished  in  their  mechanical  properties.  Diversity  of  mechanical  characteristics  stems  from 
the  variability  in  nanostuctural  organization  and  from  combination  of  different  types  of  covalent 
bonds  (sp2  and  sp3)  with  Van-der-Waals  interaction  between  molecular  clusters.  The  purpose  of 
the  paper  is  to  explain  the  elastic  property  diversity  in  the  frame  of  simple  mechanical  models  of 
different  nanoscale  structures  using,  in  the  main,  data  on  elastic  properties  of  well-known  carbon 
states  - graphite  and  diamond.  Analysis  of  these  models  makes  it  possible  to  follow  variations  in 
elastic  characteristic  of  carbon  states  in  transition  from  one  form  of  nanostructural  organization 
to  another. 

The  most  of  carbon  states  are  formed  by  2D  structural  units  - plane  or  curved  graphene  sheets 
or  closed  molecular  clusters  (fullerenes,  nanofilaments,  etc)  originated  from  these  sheets.  Van- 
der-Waals  forces  or  covalent  sp3  bonds  provide  3D  linkage  between  the  2D  graphite-like 
nanostructural  elements.  Diversity  of  carbon  states  actually  results  from  two  specific  peculiarities 
of  hybridized  C-C  bonds.  The  first  is  high  lability  - interatomic  distances  and  valent  angles  are 
liable  to  vary  over  relatively  wide  range.  Stressed  sp 2 bonds  cause  bending  graphite  sheets  - 
curved  atomic  sheets  are  inherent  in  many  carbon  states:  carbon  fibers,  fullerenes,  nanotubes  and 
so  on.  The  other  feature  of  carbon  is  concurrent  participation  of  different  types  of  C-C  bonds  in 
formation  of  carbon  states.  Normally  the  spatial  structure  of  carbon  materials  is  formed  by  Van- 
der-Waals  interaction  between  graphene  sheets  or  molecular  clusters.  This  interaction  results  in 
formation  of  soft  and  fragile  materials  (crystalline  graphite  and  graphite  materials,  pristine 
fullerites,  carbon  nanotube  paper).  Under  extreme  conditions  sp3  bonds  link  molecular  clusters 
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forming  hard  and  ultrahard  HPHT  fullerene  phases,  neutron  irradiated  graphite,  diamond-like 
films,  etc. 

Our  mechanical  models  for  evaluating  elastic  properties  of  solid  carbon  states  include:  ( 1 ) 
graphene  sheets  and  2D  closed  molecular  clusters  (fullerenes  and  nanofilaments)  as  elements  of 
spatial  construction,  (2)  Van-der-Waals  interaction  and  (3)  diamond-like  sp 3 bonds  as  factors  to 
assemble  these  structural  elements  into  an  integral  construction.  It  allows  discussing  the 
mechanical  properties  in  terms  of  2D  graphene  sheet  elasticity,  stiffness  of  individual  diamond- 
like bond  and  efficiency  of  Van-der-Waals  force. 

Carbon  molecular  clusters  involve  a significant  number  of  atoms  to  treat  them  as  macroscopic 
systems  and  to  speak  about  their  elastic  properties.  The  surface  elastic  moduli  CM  and  Cl2 
establish  interrelation  between  the  tensions  Tj,  ft  and  strains  £\,  £2  in  the  sheet  plane: 

Zi  = C„  £\  + C,,  £2,  T2  = C,,  -£|  + C„  £1  (1) 

Following  to  Ruoffs  [1],  the  2D  elasticity  of  graphene  sheets  is  evaluated  from  the  in-plane 
elastic  moduli  Cn  and  C12  of  crystalline  graphite  referred  to  a single  sheet.  The  sheet  density  N 
in  a graphite  stack  is  derived  from  the  interlayer  distance  do.  Using  experimental  data  (Cn 
=1020GPa,  C|2=  180  GPa  and  do  = 0.335  nm  [2])  we  arrive  at  the  surface  moduli  of  a graphene 
sheet: 


C,,=  Cn  -do-  341,7  Pam,  c,,=  Cn-do  = 60,3  Pa-m  (2) 

The  moduli  Cn  and  C,,have  been  used  in  [1]  to  estimate  the  bulk  elastic  modulus  fic«>  of  a single 
C«>  molecule: 


Bcm  = f +<? J = f (C, , + Cvy  A.  - 800  GPa , (3) 

/?C6o  = 0,35  nm  is  the  fullerene  radius.  The  off-plane  elastic  characteristics  (bending  resistance) 
of  the  graphene  sheet  could  not  be  estimated  from  properties  of  graphite.  Response  of  an  atomic 
layer  to  bending  could  not  be  considered  as  ordinary  reaction  of  a macroscopic  body: 
compression  of  the  area  below  the  neutral  line  and  dilatation  of  the  upper  layers.  To  evaluate  the 
response  theoretical  computation  has  been  performed  in  the  frame  of  molecular  dynamics 
methods  [3,4]  for  spherical  configuration  of  a graphene  sheet  because  of  our  special  interest  to 
elastic  properties  of  superhard  HPHT  phases  of  C60.  The  bending  stiffness  k of  a graphene  shell 
has  been  introduced  to  describe  correlation  between  the  point  force / and  the  displacement  ARcso 
of  the  shell  induced  this  point  load: 

f-kARCm/RCm.  (4) 

The  calculated  value  of  k is  57. 6 1 0'8  N. 

Elasticity  of  a single  sp 3 bond  is  characterized  by  stiffness  Ksn>  of  elementary  tetrahedron 
associated  with  the  bond  with  respect  to  the  force/acting  along  the  tetrahedron  axis: 

f = Kip  £„=Ah/h  (5) 
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The  applied  force  induces  the  change  Ah  of  the  tetrahedron  height  h that  we  consider  as  the 
length  of  the  intercluster  link.  Value  of  K^>  may  be  evaluated  from  magnitude  of  the  elastic 
modulus  Cf 1 > of  diamond  along  the  <1 1 1>  direction.  Diamond  lattice  is  treated  as  an 
arrangement  of  parallel  layers  consisted  of  diamond  tetrahedrons  oriented  along  <1 1 1 > direction. 
Normal  stress  applied  to  the  { 1 11 } plane  provides  the  force/=  er-  So  acting  on  an  individual 
tetrahedron;  So  = V3 / 4 • a1  is  area  of  the  {111}  plane  associated  with  a single  tetrahedron,  a is 
length  of fee  diamond  lattice.  The  deformation  £h  of  the  elementary  tetrahedron  is  equal  to 
deformation  induced  by  the  applied  stress  in  diamond 

eh=e  = 07 Ct  ■ (6) 

Starting  from  the  value  of  C^ul> = 1 200  GPa  [5]  the  single  sp3  bond  stiffness  has  been  calculated: 

Kv>  = f/eh=Cflt>-S„=^-ai-C^i'>=  6.58-1  O'*  N.  (7) 

Efficiency  of  Van-der-Waals  interaction  is  taken  from  the  graphite  elastic  modulus  C33 
responsible  for  interlayer  interaction  between  atomic  sheets  in  crystalline  graphite. 

Elementary  structural  units  described  above  are  employed  for  the  assessment  of  elastic 
properties  in  diverse  carbon  states.  The  attention  is  focused  on  evaluation  of  the  bulk  elastic 
modulus  B (inverse  compressibility).  The  bulk  modulus  is  a natural  elastic  characteristic  of  a 
matter.  The  bulk  modulus  can  be  evaluated  in  the  frame  of  the  models  for  all  feasible  carbon 
states.  Estimation  of  shear  elastic  characteristics  requires  more  detailed  description  of 
interatomic  interaction. 

Let  start  the  consideration  from  soft  and  fragile  carbon  states  that  are  formed  by  plane  or 
curved  graphene  sheets  linked  by  Van-der-Waals  forces.  Crystalline  and  turbostratic  graphite, 
pristine  fullerite  and  nanotube  materials  represent  this  kind  of  the  structural  organization.  Weak 
molecular  forces  between  sheets,  graphite  clusters  or  fullerene  molecules  give  the  main 
contribution  into  elastic  response  of  a matter.  In  crystalline  graphite  Van-der-Waals  interaction 
hold  together  the  stack  of  two-dimensional  graphene  sheets.  Hydrostatic  pressure  p applied  to  a 
graphite  sample  above  all  else  induces  change  of  the  distance  between  graphene  sheets.  So  the 
bulk  modulus  B is  close  to  the  interlayer  component  C33  of  the  elastic  matrice  of  graphite: 

5gr « C33  = 36  GPa. 

In  turbostratic  graphite  and  molecular  crystals  of  pristine  fullerite  weak  Van-der-Waals  forces 
are  responsible  for  interaction  between  graphite  clusters  or  fullerene  molecules  in  all  directions 
uniformly.  The  bulk  modulus  B in  these  materials  should  be  approximately  3 times  smaller 
comparing  with  the  value  of  B in  crystalline  graphite  because  hydrostatic  pressure  is  able  to 
change  intercluster  distances  in  three  orthotropic  directions.  For  a cubic  molecular  crystal  of 
pristine  fullerite  the  Hook  equation  gives  values  of  strains  £ (i  =1 ,2,3)  induced  by  hydrostatic 
pressure  p:  £.~  p / C33;  here  the  unidirectional  elastic  modulus  of  fullerite  is  approximated  by  the 
interlayer  module  C33  of  graphite.  The  bulk  deformation  £V  determines  the  value  of  the  bulk 
modulus  of  crystalline  Cm’- 

£v  = AVIV  =’£\  +£2  +£2  = 3 -p  / Cy\\  Bprjs  = C33/3  = 12  GPa.  (8) 
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The  theoretically  estimated  value  of  B is  in  good  agreement  with  the  experimentally  obtained 
magnitude  of  1 1 GPa.  The  little  difference  caused  by  curvature  of  fullercnc  molecules  because 
the  basis  of  our  estimations  is  the  molecular  interaction  between  the  plane  graphene  sheets. 

The  next  group  of  carbon  states  includes  substances  in  which  covalent  sp 3 partly  replace  the 
Van-der  Waals  interaction  between  molecular  clusters.  The  substitution  results  in  progressive 
stiffening  carbon  materials  as  increasing  quantity  of  up3  covalent  bonds.  Mixture  of  sp2  and  sp'  is 
characteristic  for  many  graphite  and  diamond-like  materials.  But  our  consideration  will  be 
focused  on  high-pressure  high-temperature  (HPHT)  phases  of  polymerized  fullerite.  Under 
HPHT  conditions  fullerites  Cco  and  C70  undergo  an  intricate  series  of  successive  transformations 
caused  by  covalent  sp3  bonding  of  fullerene  clusters.  Mechanical  properties  of  HPHT  phases 
follow  structural  transformations.  Polymerized  fullercnes  demonstrate  a wide  spectrum  of  phases 
with  diverse  nanostructure,  elastic  properties  of  that  are  available  [6]  to  compare  with  our 
modeling  estimations. 

At  first  stages  of  polymerization  fullerene  balls  are  linked  together  into  linear  chains  or  planar 
structures  [7]  by  pairs  of  .sp3  bonds  (2+2  cycloaddition  mechanism).  Molecular  chains  and  planes 
are  held  together  by  means  of  Van-dcr-Waals  forces.  Chains  of  polymerized  Cm  arc  aligned 
parallel  into  orthorhombic  crystalline  structure,  planes  provides  tetrahedral  or  rhombohedral 
structure  of  2D  polymerized  Cfin.  Van-der- Waals  interaction  between  chains  or  planes  of 
polymerized  balls  is  much  weaker  than  sp2  and  sp3  bonds  that  form  individual  chains  of  fullerene 
molecules.  Compressibility  of  the  ID  or  2D  polymerized  states  is  determined  by  changing 
distances  between  chains  and  planes  in  directions  perpendicular  to  the  direction  of  alignment. 
Their  bulk  elastic  moduli  B should  be  close  to  the  value  specific  for  pristine  crystalline  fullerite 
or  graphite  (10  + 40  GPa).  The  real  values  of  B depend  on  topology  of  material  structure.  In 
pristine  Cm  Van-der-Waals  interaction  is  responsible  for  all  of  three  deformation  ft,  £2  and  ft 
induced  by  hydrostatic  pressure  along  the  coordinate  axes.  Elastic  response  along  each  of  the 
axes  is  described  by  the  elastic  coefficient  C,  which  is  close  to  the  interlayer  elastic  parameter 
C33  of  graphite:  C ~Cn.  Because  the  bulk  deformation  £v  involves  all  three  values  £j,  e2  and  ft: 

£ v = *ft  +ft  +ft;  £\  ~ ft  — ft  ~ p / C, 

we  have  arrived  at  the  value  B,,„\,  = C/3  = 12  GPa  in  (8). 

In  ID  polymerized  Cr.n  hydrostatic  pressure  produces  two  strain  ft  and  £2  only,  the  third  strain 
ft  along  the  direction  of  molecule  alignment  vanishes  (£3  « £1.2).  So  the  bulk  elastic  modulus 
Bid  of  ID  polymerized  Cm  is  by  1,5  time  higher  than  the  modulus  of  pristine  Cm: 

B,d  = pl(£  1 +£2)  = C/2  = (3/2)  B,lni„  - 17  GPa  (9) 

In  2D  polymerized  state  nanostructure  is  easily  deformed  only  along  one  direction  - 
perpendicular  to  the  planes  of  covalently  bounded  fullerene  molecules  (£2.3  « ft).  The  sp 3 
covalent  bonding  makes  2D  polymerized  Cm  3 times  stiffer  comparing  with  pristine  fullerene: 

B2 d = p/ft  = C = 3 Bprist-  ~ 35  GPa  ( 1 0) 

Elastic  modules  of  ID  and  2D  polymerized  Cm  have  been  measured  by  the  pulse  time-resolved 
method  of  acoustic  microscopy  [6],  The  experimental  values  of  the  bulk  elastic  modulus  B for 
ID  and  2D  polymerized  Cm  fullerites  (Bid  = 27  GPa  and  B2d  = 45  GPa)  are  close  to  our 
estimations  - 17  and  35  GPa  respectively.  Divergence  of  the  experimental  and  calculated  values 
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shows  our  calculations  should  be  considered,  as  a qualitative  way  for  estimating  because  real 
features  of  interaction  between  polymerized  clusters  has  not  been  taken  into  account. 

Principally  new  situation  arises  when  sp'  covalent  bonds  pack  up  fullerene  molecules  into  a 3D 
system.  This  type  of  Cm  polymerization  is  realized  at  quite  high  pressure  (9  + 12  GPa)  and 
temperature  (higher  1000  °C).  Elastic  properties  of  3D  polymerized  fullerites  have  been 
measured  by  acoustic  microscopy  technique  [6]  and  Brillouin  scattering  method;  it  was  been 
shown  the  hardness  and  bulk  modulus  of  some  3D  polymerized  phases  are  compared  with  these 
parameters  of  diamond.  It  necessary  to  underline  3D  polymerized  phases  in  practice  are  observed 
in  crystalline  (actually .polycrystalline)  or  partly  disordered  states,  A few  theoretical  models  of 
3D  crystalline  structure  of  Cm,  fullerites  has  been  proposed  [8,9], 

For  our  estimations  we  employ  the  simplest  model  of  a cubic  3D  structure  formed  by  fullerene 
molecules  that  are  held  together  with  double  sp 3 covalent  bonding  (2+2  cycloaddition)  in  three 
peipendieular  coordinate  directions  jq,  X2,  x3.  The  primary  unit  of  this  3D  structure  is  a fullerene 
molecule  plus  six  perpendicular  double  sp3  bonds.  Elastic  response  of  the  structure  to  a stress 
applied  along  any  of  the  coordinate  axes  is  determined  by  the  combined  stiffness  K of  the  series 
system  of  elastic  elements  involved  a pair  of  covalent  sp3  bonds  and  bending  resistance  of  the 
elastic  sphere  to  a point  force.  Elastic  characteristics  of  these  elements  have  been  described 
above.  Using  the  bending  stiffness  & of  a graphene  shell  and  stiffness  K,p  of  a single  sp 3 bond  the 
expression  for  combined  stiffness  of  the  primary  structural  unit  is  written: 


K = 


^-•(2 K ,)+!*-■  k ’ 

I ' sp  * / 1 


(11) 


U - fullerene  ball  diameter  (~  ik  [10]),  h - sp 3 tetrahedron  height  (~  2A),  L - lattice  parameter 
of  3D  polymerized  C«)  (~  8,5  A).  According  to  our  previous  estimations:  k = 57.6- 10'*  N and  KSI? 
= 6.6  10'8  N.  The  resulting  value  of  the  combined  stiffness  K=  31.3  N gives  an  elastic 
characteristics  of  the  3D  polymer  Cm  structural  unit.  To  pass  to  material  properties  we  should 
transfer  from  the  stiffness  K to  the  elastic  modulus  Cu,  characterizing  one-dimensional  elastic 
response  of  a matter  dividing  the  stiffness  K by  area  S3d  = L 2 of  the  lattice  unit: 

Cu  = K/S}d  = 450  GPa.  (12) 

Estimation  of  the  bulk  elastic  modulus  B3d  for  the  simple  cubic  3D  structure  (ft-ft-Q- p/C n)  is: 

B3 d = p/£V  = Cjj/3  = 130  GPa  (13) 

The  value  of  £3 D has  been  calculated  in  approximation  Ci  1 » C12.  It  gives  the  minimum  of  the 
value.  Approximation  Cu  = Ci2  leads  to  the  maximum  value  of  Bid  ~ 400  GPa.  For  the  most 
solids  Cu  ~ 2 tn  and  S3D  is  about  260  GPa.  Experimental  values  of  bulk  elastic  modulus  B3d  ~ 

1 60  + 450  GPa  [6]  are  in  a good  agreement  with  our  calculations. 

Experimental  data  point  to  existence  of  HPHT  fullerene  states  with  values  of  elastic  moduli 
that  are  much  bigger  values  calculated  above.  To  explain  these  experimental  results  we  propose  a 
new  model  of  structural  organization  and  mechanical  properties  of  HPHT  phases.  At  very  high 
pressure  (1  l+13GPa)  and  temperature  (1200+2000  °C)  fullerenes  are  linked  through 
incorporation  of  sub-units  of  one  cluster  into  cages  of  others  (reaction  of  polycondensation  [11]), 
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The  cluster  coalescence  results  in  formation  of  rigid  and  tight-packed  3D  graphene  sheet 
network.  We  model  the  structure  as  a 3D  system  of  mutually  intersected  graphene  sheets  packed 
parallel  to  x-,  y-,  z-directions  and  connected  one  to  another  by  common  atoms  along  lines  of  their 
intersections.  The  model  seems  to  be  a good  approximation  to  estimate  elastic  response  of  ultra- 
hard  fullerite  materials. 

The  elastic  response  of  the  model  system  is  expressed  through  the  surface  elastic  moduli  c„ 
and  cp  of  a single  graphene  sheet  and  density  N=  1 //>  of  graphene  sheet  package  along  each  of 
the  principal  directions  of  packing  (/>-  is  a mean  distance  between  neighbor  atomic  sheets  along 
a given  direction).  Hydrostatic  pressure  p induces  principal  strains  ft.  ft  and  ft  along  the 
principal  axes.  In  each  of  three  orthogonal  directions  the  elastic  response  is  caused  by  in-plane 
deformation  of  two  orthogonal  systems  of  intersecting  sheets.  For  direction  1 : 

p = 2-Cu-l/hel+Ct:\/b(el+£,);  (13) 

the  same  equations  may  be  written  for  direction  2 and  3.  This  system  of  equations  makes  it 
possible  to  express  the  bulk  deformation  e v through  the  pressure  p and  to  obtain  the  equation  for 
the  bulk  elastic  modulus  B of  the  structure: 


B = 2/3-(Cn  + Cn)-{d/b),  (14) 

where  C\\,  Cn  - in-plane  elastic  moduli  of  graphite;  d and  b - interlayer  distances  in  graphite 
and  in  the  model  system  accordingly.  When  the  average  size  of  the  cluster  units  equals  to  the 
fullerene  radius  Ru  = 0.355  nm  the  expression  (14)  coincides  with  Rouffs’  calculations  for  the 
bulk  elastic  modulus  of  a free  fullerene  molecule:  B = 800  GPa.  The  ultimate  value  of  B = 

950  GPa  arises  when  the  average  interlayer  distance  b coincides  with  the  one  in  graphite. 
Calculated  values  are  able  to  explain  extremely  high  values  of  the  bulk  elastic  modulus  of 
ultrahard  fullerites  B = 550  + 800  GPa  having  been  observed  in  experiments  16]. 
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